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We investigate the polarization modulation properties of a variable-delay 
polarization modulator (VPM). The VPM modulates polarization via a 
variable separation between a polarizing grid and a parallel mirror. We find 
that in the limit where the wavelength is much larger than the diameter 
of the metal wires that comprise the grid, the phase delay derived from 
the geometric separation between the mirror and the grid is sufficient to 
characterize the device. However, outside of this range, additional parameters 
describing the polarizing grid geometry must be included to fully characterize 
the modulator response. In this paper, we report test results of a VPM at 
wavelengths of 350 /zm and 3 mm. Electromagnetic simulations of wire grid 
polarizers were performed and are summarized using a simple circuit model 
that incorporates the loss and polarization properties of the device. © 2011 
Optical Society of America 

OCIS codes: 350.1270, 120.5410, 230.4110, 240.5440, 050.6624. 
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1. Introduction 

Astronomical polarimetry in the far-infrared through microwave portion of the electromag- 
netic spectrum is a useful tool for probing the physics of interstellar dust, investigating the 
role of magnetic fields in star formation, and characterizing the radiation from the early uni- 
verse. In each of these cases, polarization modulation is an important element of instrument 
design because it enables a precision measurement by encoding the polarization information 
and thus separating it from the typically larger unpolarized background signal. 

Variable-delay polarization modulators (VPMs) operate through the introduction of a 
controlled, variable phase delay between two orthogonal linear polarizations. This is ac- 
complished in recent work by placing a polarizing grid in front of a mirror and varying 
the separation [1,2], although other related architectures have been used for polarization 
modulation via introduction of a variable electrical delay [3,4]. For a single frequency, the 
polarization transfer function of the VPM is a sinusoidal function of the phase delay be- 
tween the two orthogonal polarizations. The key to understanding the polarization transfer 
function is to determine how the phase delay is related to the grid-mirror separation, since 
the latter is the quantity that is directly measurable in an instrumental setup. In practice, 
the polarization response of a specific modulator can be measured to produce a template 
for use in reducing astrophysical data. As such, the work presented here is not required for 
demodulating a time-ordered data set. The utility of the model described here is that it al- 
lows calculation of the achievable optical response. It provides both guidance regarding the 
modulators sensitivity to design parameters and validation metrics for the observed instru- 
mental performance. In addition to VPMs, grid-mirror systems have been implemented in a 
variety of other applications [5-13]. Elements of the analysis presented here are applicable 
to these systems as well. 

In this paper, the polarization transfer function for a single modulator from the 
Hertz/ VPM instrument is measured at 350 /im and 3 mm. The observed response is well- 
characterized by a transmission line model. This work represents a refinement of previous 
results [2] in which it was qualitatively shown that VPMs modulate polarization with rea- 
sonable efficiency and without introduction of major artifacts. Here, we present new lab tests 
of Hertz/ VPM collected after improving the optical alignment. We demonstrate quantita- 
tive agreement between the measured polarization transfer function and a transmission line 
model. Also, we show that the overall modulation efficiency of Hertz- VPM is within ~2% 
of the expected value at 350 fim and in closer agreement at 3 mm. The organization of the 
paper is as follows: In section 2, we briefly review the principle of operation of the VPM. In 
section 3 we describe the circuit model for the VPMs. In section 4, we discuss the current 
state-of-the-art for analytical models used to calculate the circuit parameters and include 
numerical simulations to overcome their limitations. In sections 5 and 6, this model is applied 
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to measurements of the VPM, and the results are discussed in section 7. 



2. The Variable-delay Polarization Modulator (VPM) 



Input Port Output Port 



Fig. 1. The VPM consists of a polarizing grid placed in front of and parallel to 
a planar mirror. The polarization parallel to the grid wires is reflected by the 
grid. The orthogonal linear polarization passes through the grid and is reflected 
off the mirror. The two components are recombined at the output port with 
a relative phase delay that is dependent upon the grid-mirror separation, d. 
The wire grid spacing or pitch is indicated by g. 



The VPM consists of a polarizing grid positioned parallel to and in front of a mirror. A 
change in grid-mirror separation corresponds to a change in introduced phase between two 
linear orthogonal polarizations. The device is show in Figure 1. Choosing coordinates such 
that Stokes Q gives the difference between the polarization states parallel and perpendicular 
to the VPM wires, the polarization transfer function can be expressed as 

U' = Ucosd + Vsin6. (1) 

Here U and U' are the input and output Stokes U parameter, and V is the input circular 
polarization. Here 5 is the electrical phase delay between the polarized component transmit- 
ted by and that reflected by the grid upon recombination at the output port of the device. 
This is the phase of interest when using the VPM as a modulator. 

In the limit in which the wavelength is much larger than the length scales that characterize 
the local grid geometry, the VPM phase is proportional to the path difference between the 
two polarizations, 

£oo = — r- cos6>, (2) 
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where 9 is the incident angle, as shown in Figure 1. Here, the oo subscript indicates that 
this relationship holds in the long wavelength limit. In many polarizing grid applications the 
size of the wires relative to the wavelength is large enough to produce measurable deviations 
from the phase relation given in Equation 2. 

We assume throughout this paper that the grid wires are uniform, infinite in extent, and lie 
in a plane. Diffraction grating lobes are negligible when the wire radius, a, and grid constant 
or center-to-center wire pitch, g, are much smaller than the wavelength of the incident 
radiation. In this limit, the structure has a homogeneous response for each polarization. 

3. A VPM Model 

At the highest level, our goal is to describe the polarization modulation properties of the 
VPM. To do this, several details are worth briefly considering before we proceed. Both Jones 
matrices [14, 15] and Mueller matrices allow one to propagate the polarization properties 
through a system having no reflections between elements. These two descriptions can be 
shown to be mathematically equivalent [16]. A polarization modulator, such as the VPM, 
that inherently relies on the interference of the fields in the device is not directly amenable 
to this analysis approach. To incorporate the influence of multiple reflections in such a 
system, the U ABCD" matrix approach can be employed. The ABCD matrix is also known 
as a transfer, transmission line, chain, or characteristic matrix [17-20] in the microwave 
and optical literature. Physically, this matrix formulation links the propagation of the input 
and output electric and magnetic fields through a given component. Once the modulator 
is modeled with these tools, it is straightforward to calculate a Jones matrix for the entire 
system [21], assuming that the input and output ports are matched. 

From an instrumenation perspective, an "ideal" polarizing grid perfectly separates or- 
thogonal linear polarization components with a response that is independent of frequency. 
That is, the grid perfectly transmits the polarization component having the electric field 
perpendicular to the grid wires and perfectly reflects the component having the electric field 
parallel to the grid wires. The former corresponds to zero impedance contrast between free 
space and the grid; the latter corresponds to an infinite impedance contrast. For physically 
realizable grids, the achievable impedance contrast for each of the two polarizations de- 
pends on the electromagnetic properties and geometry of the structure. In the limit the wire 
radius and separation are finite compared to the wavelength of the incident radiation, the 
response becomes dependent on frequency. Motivated by the observation that the wire grid is 
a polarization- and frequency-selective surface, we model the structure using a transmission 
line approach. 

We model the polarizing grid using a circuit representation [22] of the VPM as shown in 
Figure 2. We normalize all impedances to that of free space, Z = 377 Q per square. For the 
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polarizing grid, the transmission and reflection of the electric field component parallel to the 
grid wires can be represented by an inductive circuit. The circuit reactances used are defined 
in Figure 2. We compute the ABCD matrix from these as follows: 
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(3) 



Likewise, the grid's effect on the polarization component having its electric field perpendic- 
ular to the grid wires can be represented by a capacitive circuit and is given by 
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Fig. 2. In the limit the wavelength is large compared to the wire pitch, the VPM 
can be modeled by two independent circuits. An inductive circuit is used for 
the polarization component having the electric field parallel to the grid wires 
(top). A capacitive circuit is used to model the polarization component having 
the electric field perpendicular to the grid wires (bottom). 



Each of these matrices can be cascaded with the transmission matrices for a free space 
delay and a short. The resulting matrix is employed to calculate the reflection coefficients 
for each polarization, T» and r_|_, as a function of frequency. The Jones matrix for the VPM 
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Fig. 3. The geometry for the general polarizing grid lying in the x-y plane is 
shown. 




can then be expressed as 

7 -(?r"l- 

The density matrix describing the polarization state may be written as a linear combination 
of Stokes parameters (I, Q, U, V) [15]. 

n\ / n 1 \ / n „■ \ 

(6) 

The density matrix of the incoming light can be mapped to that of the outgoing light, D', by 
using the system's Jones matrix D = j' D J. This model is applicable for normal incidence. 
In the more general case of non-normal incidence, it is convenient to parameterize the model 
in terms of the incidence angle, 9, and the angle between the grid wires and the plane of 
incidence, <j). See Figure 3. This general case can be calculated using the above formalism 
along with a procedure similar to that described by Goldsmith [18] where the free-space TE 
and TM polarizations are projected onto the grid wires. The polarization separation then 
occurs in the basis of the wires, as done above. 

If the incidence angle is non zero, there are two effects. First, the length of the grid-mirror 
separation in the transmission line is multiplied by cos 8. Second, it is necessary to take into 
account the differences at the grid interface for the TM and TE polarizations. For the case 
of the transverse magnetic (TM) polarization mode, the effective impedance is altered from 
the free-space value by a factor of cos^. This is due to the fact that the electric field that 
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is projected onto the plane is smaller than the freely-propagating value. For the transverse 
electric (TE) polarization, it is the magnetic field vector that is projected into the plane and 
so the impedance differs from that of free space by a factor of 1/ cos 6. 

A simple case that one can consider for non-zero incidence angle is that where the wires 
are either parallel or perpendicular to the plane of incidence. In this case, each of the TE 
and TM modes correspond to only one of the circuits described above. For example, if 
the wires are parallel to the plane of incidence, the TE and TM modes correspond to the 
capacitive and inductive circuit elements, respectively. For a general grid rotation angle 0, 
four circuit models must be considered utilizing all of the combinations of (TM,TE) and 
(inductive, capacitive). After the transmission matrices are found for the two transmission 
lines in Figure 2, scattering parameters can be determined. Once the scattering matrices are 
found for the VPM in each of these four cases, we can combine the results into a single Jones 
matrix, 
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(7) 



This Jones matrix maps the incident electric field ^ E™ E TE J to the reflected field 

y E™' E TE ' j . The resultant polarization transformation for the system can be com- 
puted from D = J D J. 

4. Circuit Parameter Values 

The circuit topology described above is physically motivated; however, circuit parameters 
for the grid are required to analyze the response of the VPM. Many treatments found in 
the literature for wire grid polarizers consider a simple analytically-treatable limiting case 
[8,22-24]. These models assume no diffraction, A > 2g, little or no azimuthal dependence of 
current density on the wires, a < A, and that the grid filling fraction is small, 2a/ g < 1/27T. 
However, at millimeter and submillimeter wavelengths, realizable grid structures often violate 
the third condition as indicated in Figure 4. A survey of practical grid polarizers indicates 
that filing fractions in the range 0.3 < 2a/ g < 0.8 are available and achieve high polarization 
isolation. Such grids are outside of the range of applicability of the perturbation-based models 
referred to above. An interesting extension of this approach utilizing higher-order impedance 
boundary conditions applicable for sparse grids has been explored [25]. Green's function 
and lattice-sum [26-28] methods enable more accurate treatments without a need for these 
approximations. Here we employ a hybrid numerical approach that allows us to analyze the 



7 



behavior of wire grid polarizer used for this work. In this section, we first revisit the classic 
analytical treatment of polarizing grids and motivate the need for a more rigorous treatment. 



4- A. Sparse Grid Approximation 

As an example of commonly-employed analytical methods, we explicitly give the expressions 
for the circuit elements from Marcuvitz [22] ; however, this is representative of the approxima- 
tions used in quasioptical treatments of grids. In the following, all reactances are normalized 
with respect to the impedance of free space. Given a wire radius a and a separation or 
"pitch", g, the circuit elements for the polarization parallel to the wires (the "inductive" 
mode; Marcuvitz §5.21, "Inductive Posts") are given by 




The series for the inductive reactance in Equation 9 is rapidly convergent. In the limit of 
applicability, this series can be approximated by ~ 0.6 ■ (g/X) 2 and is a small perturbation. 
As a result, the factor in curly braces, which we identify with the quasi-static inductance, is 
essentially frequency- independent. 

For the polarization having the electric field perpendicular to the wires (the "capacitive" 
mode; Marcuvitz §5.13, "Capacitive Posts"), the reactive circuit elements are given for a 
symmetrically-centered capacitive post [22]: 
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The circuit parameters given above are plotted as functions of the geometric filling frac- 
tion, 2a/ g, in Figure 4. In plotting the circuit parameters, we multiply the reactances by 
a factor of X/g or g/X for inductive and capacitive components, respectively. This sepa- 
ration between wavelength-dependence and the structure's geometry is anticipated in the 
quasistatic limit (i.e. the limit in which each circuit element in the model can be treated as 
a discrete frequency-independent component). The triangles at the bottom of the figure in- 
dicate typical values of grids employed in the submillimeter and millimeter, and the vertical 
gray line marks the position at which 2a /g = l/ir. We note that the expression for the shunt 
inductance is zero at this point and negative beyond it. This is a manifestly unphysical result, 
as the wires do not actually become "capacitive" for the parallel polarization, but approach 
the inductance of a corrugated conductive plate. Thus, a different approach is required in 
this regime, as the limit of applicability indicated by Marcuvitz is violated. 
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Fig. 4. The Marcuvitz circuit parameters are plotted as functions of the ge- 
ometric filling fraction, 2a/ g. The vertical gray line is shown at 2a/ g = 1/n. 
The shaded region indicates the region in which the Marcuvitz approxima- 
tion ceases to hold [22]. Triangles indicate typical values of 2a/ g for grids 
constructed for submillimeter and millimeter use. 



4-B. Numerical Simulations of Grid Performance 

In order to understand the inductance and loss of wire grids, we perform numerical simula- 
tions as a function of filling fraction over the range 0.02 < 2a/ g < 1. The numerical grid sim- 
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ulations were carried out with CST (Computer Simulation Technology) MICROSTRIPES™ 
Time-Domain Transmission- Line Matrix (TLM) solver and Ansoft HFSS (High Frequency 
Structural Simulator). See Figure 5 for the electromagnetic configuration and the boundary 
conditions used for modeling each polarization. We consider the grid's absorption from the 
structure in the limit where the field penetration depth is small compared to the wire di- 
ameter and the incident wavelength is greater than two times the wire pitch. In this limit 
the higher order Floquet harmonics or grating lobe responses are absent, and the fields are 
quasi-static in nature and a perturbation on the lossless case. 

Wheeler's incremental inductance rule [29] can be applied to compute the resistive losses 
from the wires in the grid. This approach implicitly assumes that the current distribution 
on the conductors does not vary appreciably over distances comparable to the thickness of 
the wire and is a convenient parameterization. This condition is satisfied for a wire having 
a well-developed skin effect, 5 <C a, in the single mode limit. The anticipated form can be 
expressed as a product of two separable functions: the first dependent on the grid filling 
factor, 2a/ g, and the second a function of skin depth over the wavelength, 5/X. Guided by 
this physical insight, we proceed to model the grid in the single mode approximation to 
further explore the scaling properties of the grid inductance and loss. 



Unit Cell: 
Inductive 

2a 

d «2a<g<XI2 
Perfect-H 




Unit Cell: 
Capacitive 

2a- 

d«2a<g<XI2 
Perfect-E 




Perfect-E 



Perfect-H 



Fig. 5. The unit cells used for the simulation in both the inductive (left) and 
capacitive (right) modes are shown. The large arrow in each case indicates the 
direction of the incident electric field in the simulation. Perfect- "E" and - "H" 
indicate the use perfect electric and magnetic mirrors on the boundaries of 
the unit cell, respectively. Plane wave illumination of the unit cell with these 
boundary conditions allows representation of an infinite grid for A > 2g. 



The power reflection coefficient, R, from the impedance presented by the shunt reactance 
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of the grid, Z L = R L + iX L , is: 

R = IiSii h 



1 + 2Z T 



[i + 2R L y + (2x L y 



(14) 



where all impedances are normalized to Z . The resulting reflection phase is 



r = vr - arctan ( ^gj = " " ^ (lT^) " (15) 

For a near optimal grid geometry, the impedance contrast with respect to free space, 2Zl <C 
1, and the resultant reflection phase is approximately n. Similarly, the power transmission 
can be expressed as 



T= \S- 
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and the power absorption, A = 1—R — T = 1 — |Sii| 2 — |S*2i| 2 , is computed from conservation 
of energy. Here, Sij are the complex scattering parameters for the structure. A corresponding 
set of expressions can be written for the capacitive case. We note the ratio of the absorptance 
over the reflectance can be expressed as 

A 4R S 4tt<5 

- = AR L = ri— = v— , (17) 

where we recall from Figure 2 that Rl is the inductive mode's resistance normalized to 
the impedance of free space. The bulk surface resistivity of the metal is Rs = 1/crS, r\ is 
the grid loss efficiency, and A is the observation wavelength. The field penetration depth, 
5 = y/2/Jww^ is computed from the metal's permeability, /i, electrical conductivity, a, and 
the observation frequency, uj. It is informative to consider the dimensionless ratio that occurs 
in the circuit model above in terms of this parameterization: 



The expression on the right is related to the Hagen- Ruben emissivity formula [30]. This can 
be derived from Fresnel coefficients and identified as the emissivity from a bulk metal at 
normal incidence. This approximation is valid in the limit of low emissivity. Equation 18 
links the circuit theory to the interaction between the materials and the electromagnetic 
fields. We see that the loss efficiency, rj = (A/R)(4tt5/X)~ 1 , is the loss of the grid relative to 
a flat sheet made of the same material. We compute the scattering parameters as a function 
of grid filling factor and frequency and derive equivalent circuit parameters. In Figure 6 we 
show the grid loss factor, the grid shunt reactance, and the grid reflection phase appropriately 
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Fig. 6. (TOP) The grid loss factor, 77, is plotted as a function of geomet- 
ric filling factor for a wire grid. Both HFSS and Microstripes simulations 
are shown for the cases of the electric field perpendicular to and parallel to 
the grid wires. The filling factors for typical polarizing grids are plotted for 
comparison. (MIDDLE) The reactance for the inductive mode is shown as 
a function of the geometric filling factor. A useful interpolation function is 
X n = 0.51 - 1.191n(27ra/#) + 0.53(ln(27ra/#)) 2 + 0.11(ln(27ra/#)) 3 . For com- 
parison, the Marcuvitz model is shown as well as the measured inductance for 
the grid. (BOTTOM) The normalized reflection phase of the inductive mode 
is plotted as a function of geometric filling factor for the HFSS simulation. In 
this case, we find <fi n = 1 — A/3.545(2a/g). In each of the panels, the condition 
2a/ g = I/71", is denoted by a dashed vertical line. Triangles indicate typically 
manufactured grid geometries. 
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normalized to remove the wavelength dependence as a function of grid filling factor. This 
scaling is anticipated and observed in the quasi-static limit for these circuit elements. 

For the perpendicular illumination case, we find a grid loss factor, rj± = 2a/ g. For the 
parallel illumination case, some care is required - the conductive wires effectively define a 
finite thickness aperture. For this configuration, the influence of the waveguide cutoff has a 
non-negligible influence on the grid's absorption properties, and we find its influence results 
in a grid loss factor of 7711 = 7/tanh(77r ■ 2a/ g) for this polarization. Here, 7 = 1.15. In 
referring to Figure 6 we note that this differs from previous grid loss factors reported in the 
literature. For example, for infinitely thin strips, Ulrich [31] finds 7711 = g/2a. For small filling 
factors, this expression and the one derived above have the same slope; however, in the limit 
2a/ g — > 1, the influence of the exclusion of the fields by the grid wires has a pronounced 
influence on the response. This behavior is also observed in our numerical simulations. We 
note that within the context of the planar waveguide model [29], by symmetry, the field 
configuration for this polarization is physically equivalent to that of a microstrip via [32]. 

The numerical approach described here was pursued to overcome and understand the limi- 
tations of the perturbation expansions in aj A [22-24] and 5/ A [8]. The normalized reactance, 
X n , of the inductive component is computed in the range of interest and is plotted as a func- 
tion of 2a/ g in the middle panel of Figure 6. For the Marcuvitz model that is also shown, X n 
becomes negative for 2a / g > 1/tt which is a symptom of a failure of the asymptotic expan- 
sion and anticipated from the restrictions on the analytical expression as specified [22] . This 
generic behavior is present in other analytical treatments that are expansions in a/A [8,23,24] 
and anticipated since the derived expressions only have logarithmic accuracy. Higher order 
terms are required in this regime, e.g., see Landau §33 [33], and [34]. However, as g/\ — > 0, 
\Xl\ — > as well, and as a result, this particular modeling detail is of reduced importance. 
This analysis suggests that the observed monotonic flattening of 7711 and X n at larger values 
of 2a /g result in high performance grids; however, the reflectivity and transmissivity are 
more effectively improved by reductions in g/\. 

5. 350 /im VPM Tests 

5. A. Experimental Setup 

Preliminary characterization of the Hertz VPM instrument has been reported [2]. Our setup 
is similar to that used previously and is schematically depicted in Figure 7. A chopped black- 
body is polarized by a wire grid and focused onto the Hertz polarimeter that simultaneously 
measures horizontal and vertical polarizations. The VPM under test is placed near a pupil 
in a collimated beam within the fore optics that couple the source to the Hertz cryostat. For 
a given grid-mirror separation, the polarization is modulated by rotating a cold half-wave 
plate (HWP) in discrete steps. The aperture of the blackbody is sized to roughly fill a single 



13 



Hertz detector. Inside of the Hertz cryostat, the polarization is separated into two linear 
orthogonal components by a wire grid analyzer. Both the transmitted (T) and reflected (R) 
components are detected by bolometers. 




J Modulator Throw 

VPM 

(wires parallel to page) 

Fig. 7. The main elements of the experimental setup for the 350 /im test 
are shown [2]. Radiation from a blackbody source is polarized by a wire grid 
polarizer having wires oriented at an angle of 45° with respect to the plane 
of the page. A chopper modulates the intensity of the signal. The radiation 
is collimated prior to being processed by the VPM. Upon exiting the VPM, 
the radiation is relayed to the Hertz cryostat. Inside, the radiation passes 
through the half-wave plate and the bandpass filter before being diplexed into 
two orthogonal linear polarizations. A bolometer detects the signal in each 
polarization. 

The polarization transfer function of the VPM is measured as follows. At various grid- 
mirror separations, the output linear polarization is measured using the Hertz polarimeter. 
The normalized Stokes parameters are then reported as a function of the grid-mirror sepa- 
ration. 
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5.B. Data Reduction 

At each grid-mirror separation, the basic data analysis pipeline is similar to that described 
in [35]. However, we describe the process in some detail since there are some deviations 
that are associated with the difference between laboratory measurements and those done 
astronomically. At a single HWP position, 16 "frames" are recorded. Each frame consists of 
two co-added "chops" or on-off differences. In our laboratory setup, the chopper phase was 
observed to drift relative to the commanded TTL signal from the data system. To correct 
for this, the in-phase data and the quadrature data were both used and the signals for a 
single frame in the R and T arrays are 



ri{dj) 

ti{0Cj) 



in phase/ 

r 4 [on 



An phase/ 

t- (a,- 



quad 



tT ad (a, 



1 2 



(19) 
(20) 



Here, ctj is the value of the HWP angle at the jth HWP position. This technique mitigates 
the problems associated with phase drift in the chopper. 

The frames recorded for each HWP are averaged separately for each of the R and T signals, 
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Here, N = 16 is the number of data frames taken at each HWP position. This process is 
repeated for each of six HWP positions (nominally, 0°, 30°, 60°, 90°, 120°, and 150°). The 
relative gain between the R and T bolometers is calculated as the ratio of the sum of all the 
data in the R bolometer and that of the T bolometer: 

E 3 R(uj) 



f 



(23) 



The polarization signal is then calculated using both the R and T bolometers. This combi 
nation aids in the removal of common-mode noise in the system, 

R( aj ) - fT(a 3 ) 



S(otj) 



R(aj) + fT(aj) 



(24) 



For each grid-mirror separation, the HWP polarization signal in Equation 24 was fit to 
the modulation function 



S(a) =ei(gcos4a + u sin 4a) + ei cos la + €3 sin 2a. 



(25) 
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The parameters q = Q/I and u = U/I are the normalized linear Stokes parameters in a 
rotational basis for which q defines the difference in polarization components perpendicular to 
and parallel to the grid wires of the VPM. We have explicitly included the linear polarization 
efficiency, ex. The system was found to be dominated by drifts on the time scale associated 
with the HWP cycle. To account for this, error estimates are determined by goodness-of-fit 
of the data to Equation 25. 

For the purposes of fitting, e 2 and €3, are orthogonal to those of importance, q and u and 
thus are discarded for the primary analysis. However, to verify our understanding of the 
system, it is worth briefly exploring their physical origin that stems from well-understood 
non-ideal behavior of the HWP [36]. 

Because the zero angle of a is set by the VPM grid wires and not the fast axis of the 
HWP, e 2 and €3 are linked by a rotation to coupling to both unpolarized flux and circular 
polarization, 



Here, x is t ne angle between the VPM wires and the fast axis of the HWP. The factor e' 2 
is the coupling between unpolarized flux and the polarization signal due to bi-attenuance in 
the HWP. The factor e' 3 is proportional to the circular polarization (Stokes V) incident on 
the HWP. The response of a single-layer HWP, being a function of the birefringence of the 
crystal (An) and its thickness, is only ideal at a single frequency. Deviation from perfect 
matching between the HWP center frequency and that of the band-defining filter as well as 
the details of the input spectrum can both cause some leakage between circular and linear 
polarizations. 

The VPM provides a means for separating the two effects. As the grid-mirror separation 
is varied, e' 3 will vary with phase delay while e' 2 will be unmodulated. An unmodulated 
offset was not observed, implying an upper limit of e' 2 < 0.01. From the modulated signal, 
we find e 3 ~ 0.2. At this level, a corresponding reduction of the HWP response to linear 
polarization is expected. We calculate this effect to be ^/l — e 3 2 = 0.98, thus resulting in an 
efficiency reduction of ~ 2%. This is consistent with previous characterization of the Hertz 
polarimeter [37,38]. In the test configuration with the VPM, the HWP was exposed to much 
larger amounts of circular polarization than in an astronomical instrumentation setting, thus 
emphasizing this otherwise small effect. 

5. C. Systematic Offset Removal 

In our test configuration, there are three aluminum mirrors between the source grid and 
the VPM having incidence angles ~ 20°. From this we estimate an induced polarization of 
< 0.3% that is primarily in the ±g direction. Instrumental polarization between the VPM 
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and half-wave plate would lead to a small change in the efficiency of the system. From 
previous laboratory measurements [2] we estimate a resulting systematic uncertainty in the 
polarization efficiency of < 1%. 

Initially we observed that Stokes q exhibited a modulation as a function of grid-mirror 
separation that had a similar modulation pattern to that anticipated for circular polariza- 
tion. This effect was observed at the 5% level. By rotating the Hertz cryostat by 90° and 
remeasuring, this effect was found to switch signs, indicating some birefringence in either the 
lenses or the pressure window in the cryostat. Since the axis of such a birefringent element 
could a priori be oriented in any orientation, a similar effect may be contaminating Stokes u. 
To cancel this systematic effect the above-mentioned symmetry, the reported values for the 
VPM transfer function were evaluated by averaging each data point with an equivalent data 
point taken in a configuration where the cryostat was rotated by 90°. The resulting values 
for q and u are shown as a function of grid- mirror separation in Figure 8 A. 

5.D. Circuit Model Parameter Extraction 

We have applied the circuit model to data taken using the Hertz VPM instrument with a 
single VPM having wires nearly parallel to the plane of incidence. In this case, it is necessary 
to average over all relevant frequencies: 



where ip n {v) is the normalized instrument bandpass (approximately z/q = 353 /im, Au/u = 



The q and u data are fit to the VPM model using the downhill simplex method [39]. Ten 
fit parameters were varied in this minimization: the angle of incidence, an offset in the grid- 
mirror separation, a global polarization efficiency, a HWP phase offset, Xlc Xl, Rl, Xq, 
Xcci and the bandwidth stretch (k). The bandwidth stretch takes into account deviation 
from parallelism between the grid and mirror in the VPM. Such an effect would cause the 
bandpass to appear slightly wider than that defined by the bandpass filter. 

The best fit parameters for our circuit model are given in Table 1. Where applicable, these 
values are compared with the values estimated from the system setup and the derived circuit 
values for the grid (a = 12.5 /im and g = 67.5 /im.) The resulting models for q and u are 
plotted along with the data in Figure 6A. The reduced x 2 value for the joint (q, u) fit is 2.1. 
Upon removal of five of the (q, u) pairs where either q or u has a significant deviation from 
the model, the x 2 is reduced to 1.2. From this we conclude that the circuit model captures 
the basic features of the polarization transfer function, however, the statistical errors are 
underestimated in a subset of the measurements. Such an effect can occur as a result of 
time-variable instabilities in the system that are insufficiently mitigated by the techniques 
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described earlier in this section. Based upon the observed performance of the system a phase 
drift in the chopper wheel for a subset of the data frames is considered the most probable 
cause. 

The resistance in the inductive mode, Rl, is higher than expected from the formulation 
given above. The expected value is based on the conductivity of gold, the coatings for the 
wires. This is an indication that the gold coatings are compromised; the fitted resistance is 
consistent with the base metal (tungsten, p = 1.8 x 10 7 Q ■ m [17]). The fitted value of Xl is 
within 30% of that derived from the numerical simulation described above. We have plotted 
this value shown in the middle panel of Figure 6. The model presented here depends on the 
ability to identify a unit cell for the grid; however, variations in wire spacing violate the 
symmetry used to calculate the properties of the grid using this technique. Examination of 
the grid using a microscope confirmed this concern. 

In Figure 8B, the features in Stokes q at grid-mirror separations of approximately 200 
and 400 /im are due to trapped resonances for the inductive mode. We note in passing that 
these features can be significantly reduced by modifying the VPM configuration. In place 
of the mirror, one could incorporate a second grid with wires oriented perpendicular to the 
first followed by beam dump. Ideally, the brightness temperature of the beam dump should 
match that of the source over the beam on the sky. This architectural change would reduce 
common-mode systematic effects between orthogonally-polarized detectors. 



Parameter 


Estimate 


Fit 


Units 


Angle of Incidence 


20 


22.83 


[degrees] 


Grid-mirror Offset 





-2.01 


[/im] 


Polarization Efficiency 


100 


98.4 


[%} 


HWP Phase Offset 





0.0125 


[radians] 


H 


1 


1.08 




Rl 


9.9 x 10~ 4t 


3.4 x 10~ 3 




x L 


0.031+ 


0.040 




Xlc 


-0.26* 


-0.27 




X c 




-0.43 




Xcc 




-265 





Table 1. Fit results for the 350 /im measurements. (^Frorn §4, *From Marcuvitz) 
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Fig. 8. (A) The q and u values were measured as a function of grid-mirror 
separation at 350 /im. (B) Measurements of the VPM at 3 mm are also shown. 
The transmission line model is in close agreement with the infinite wavelength 
approximation at this frequency; however, the difference between the two is 
measurable. 
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5. E. VPM Modulation Efficiency 

The modulation efficiency of the Hertz polarimeter was experimentally observed to be ~ 95%. 
This measurement was made with the source grid in place but in the absence of the VPM 
polarizer. This value is consistent with previous measurements of the efficiency of the Hertz 
polarimeter [2,38]. For the data reported in Table 1 and Figure 8, the reported efficiency, 
98.4%, is that of the VPM alone. The lower modulation efficiency previously reported, 85% 
[2], was found to be reproducible and induced by beam misalignments at the location of 
the VPM. The estimate reported here for the modulation efficiency of the VPM is a lower 
limit and could arise from one of several effects. First, the circuit model assumes uniform 
wire spacing and size. As such, the model represents the average circuit parameters over the 
illuminated part of the grid and does not account for variations in parasitic coupling between 
the wires that are important at high frequencies [40,41]. Treatment of this effect would require 
modifications to the simple circuit model considered here. Alternatively, and more likely, the 
non-unity efficiency factor could be an artifact of the system such as a residual misalignment 
of the beam on the VPM or residual uncertainty in the Hertz efficiency calibration. 

6. 3 mm VPM Tests 

We have also measured the VPM at a wavelength of 3 mm (g/X = 0.022.) We measured 
the response from 75-110 GHz; however, we used the 91-100 GHz portion of the measured 
spectrum to synthesize a 10% bandwidth to match the 350 /im measurement. The setup for 
this test is shown in Figure 9. An HP 8510 vector network analyzer was used to measure 
the polarization state of the radiation as it is processed by the VPM. Vertically-polarized 
radiation is emitted from port 1 and launched from waveguide to free space via a feedhorn. 
It is then reflected by a polarizing grid having wires oriented in the vertical direction. This 
grid helped to further define the polarization state of the radiation and limits the formation 
of trapped reflection resonances in the system. The unused polarization is terminated by a 
free-space beam dump. An ellipsoidal mirror is used to form a beam waist on the VPM. 
This ensures that the phase front of the radiation is approximately flat at the polarization 
modulation point. Just as for the 350 fim setup, the VPM wires were oriented such that their 
projected angle is 45° with respect to the incoming polarization state. A second, identical 
ellipsoidal mirror refocuses the beam on a second feedhorn that is identical to the first. A 
folding flat is used in this final free-space path to symmetrize the optics. The orthomode 
transducer (OMT) [42, 43] ports corresponding to the unused polarizations are terminated 
by precision waveguide loads. This mitigates the influence of polarization mode conversion 
that can result in trapped resonances in the beam waveguide structure. 

As the grid-mirror separation is varied, the VPM modulates the incoming (vertical) signal 
that we define as Stokes U in order to maintain consistency with the 350 /im data. The 
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modulated power is detected by port 2 which is sensitive to vertically-polarized radiation. 
The 800-point spectrum of S12 over the entirety of W-band (75 to 115 GHz) is recorded at 
every grid-mirror separation. Once this is finished, a 90° twist is added to the waveguide in 
port 2 to switch the port's sensitivity from vertical to horizontal polarization. The grid-mirror 
separation of the VPM is then adjusted in an identical way as previously described. 




(wires at 45 deg.) 



Fig. 9. The setup for the 3 mm VPM transmission test is shown. Radiation is 
emitted from port 1 of the HP 8510 with a vertical polarization. The radiation 
is then reflected off of the grid and collimated by an ellipsoidal mirror. At this 
point, the radiation is reflected off of the VPM. An identical ellipsoidal mirror 
follows, and a folding flat directs the beam into a second feed that is attached 
to port 2 of the HP 8510. A 90° twist is added to the waveguide to change 
the sensitivity from the vertically-polarized state to the horizontally-polarized 
state. 



6. A. Data Analaysis 

For comparison with the 350 /jm data, we wish to determine the response of the VPM over 
a 10% bandwidth fractional bandwidth. To do this, we calculate the integrated signals H(z) 
and V(z): 
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S(d)= E \ S ^)\ 2 (28) 

Au 

V(d)=Y,\ S u( d ^\ 2 - ( 29 ) 

Au 

Here, Si 2 {d, v) and S%(d, are the relevant scattering parameters for the case of horizontal 
and vertical sensitivity for port 2, respectively. 

The next step is to calculate the relative gain of the two channels. This can be done by 
noting that the single-frequency response for the two quantities are, 

H(d) = ~ + U H cos <j>(d) + V H sin 0(d)] (30) 

V(d) = ~[I V ~ U v cos<f){d) - V v sin 0(d)]. (31) 

Here, Ih/v, U H / V , and V H / V are Stokes parameters. 

For each of H(d) and V(d), we perform a linear least-squares fit. Since the single-frequency 
will not fit our broadband response precisely, we do not use this to directly calculate Stokes 
parameters. Rather, we calculate a relative gain based on the unmodulated component 
(Stokes I) from each fit, 

/ = T-. (32) 

From this, Stokes u can be calculated: 

u(d) - H{d) ~ fV{d) (33) 

u[d) ~ H(d) + fv(dy [66) 

The results of this analysis are shown in Figure 8B. In addition, we have plotted a VPM 
model using the parameters found from the 350 /xm fit, appropriately scaled by wavelength. 
We found that to fit the data required an incident angle of 19° and an offset in the grid-mirror 
separation of -30 /im. No other parameters in the model have been adjusted. For comparison, 
we have plotted the normalized Stokes parameter, u, in the infinite-wavelength limit. In this 
limit, the impedance contrast is high for the inductive mode and the dependence on the 
circuit parameters is minimal. 

7. Discussion 

Figure 10 summarizes the VPM transfer function for a single frequency for each of the 
models in this work as well as the geometric limit. To emphasize the differences in shape 
of the curves, for each case we have shifted the curve to align the first minimum to lie at 
zero delay. For A ^> a, the grid circuit is firmly in the single-mode limit. In this case, for 
plane wave illumination, Equation 2 is a good approximation for the VPM phase delay. As 
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the wire diameter becomes a finite fraction of a wavelength, the impedance contrast of the 
inductive circuit is reduced. The polarization response remains a sinusoidal function of the 
phase delay; however, the VPM reflection phase is dependent upon the details of the grid 
geometry and impedance contrast. 

8. Summary 

We have applied a transmission line model to the VPM used for the Hertz/ VPM submil- 
limeter polarimeter and have obtained estimates for key circuit parameters using numerical 
simulations. The simulations represent an improvement over analytical approximations in 
estimating the inductance of wire grid polarizers having 27ra ~ g, such as those used in this 
work. 

At 350 fim the VPM grids have wires that are ~0.07A in diameter, corresponding to a 
free-space phase shift across the width of a wire that exceeds 10°. We have shown that in 
this regime, it is not sufficient to identify the VPM phase delay with that derived from 
the path difference induced by the grid-mirror separation. The VPM transfer function is 
modified from the infinite wavelength case, both in shape and in position. This is due to a 
combination of the reactances imposed by the thickness and spacing of the wires and the 
interaction of these reactance with the variable position of the mirror. At 3 mm, where the 
impedance contrast for the inductive mode is higher, differences from the geometric phase 
delay (Eq. 2) are reduced, but evident. 
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Fig. 10. The modeled VPM reflection phase delay is shown for single- 
frequencies at 350 and 3000 /im for the grid models obtained above. In the 
limit g/X <C 1, a sinusoidal form for Stokes U is observed. As this condition is 
relaxed, the VPM reflection phase delay differs from the free-space grid-mirror 
delay. The legend on the right shows the parameters corresponding to each of 
the curves. 
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